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Abstract. Quantum Spin-Hall systems arc topological insulators displaying 
dissipationless spin currents flowing at the edges of the samples. In 
contradistinction to the Quantum Hall systems where the charge conductance 
of the edge modes is quantized, the spin conductance is not and it remained an 
open problem to find the observable whose edge current is quantized. In this 
paper, we define a particular observable and the edge current corresponding to 
this observable. We show that this current is quantized and that the quantization 
is given by the index of a certain Predholm operator. This provides a new 
topological invariant that is shown to take the generic values and 2, in line with 
the Z2 topological classification of time-reversal invariant systems. The result 
gives an effective tool for the investigation of the edge structure in Quantum 
Spin-Hall systems. Based on a reasonable assumption, we also show that the 
edge conducting channels are not destroyed by a random edge. 
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1. Introduction 

A new class of insulators has been recently found [U [21 [3l H] to possess a dissipationless 
Quantum Spin-Hall effect. Describing the structure of the edge modes in these 
systems remains an interesting issue for both fundamental understanding and potential 
applications of the Quantum Spin-Hall effect [5j [TJ [HI IHJ EHl Ell HH US] • It was argued 
in the literature that the initial Z2 topological classification proposed in Ref. [1] can 
be further refined to meet this purpuse. Not long ago, Ref. [6] introduced a new 
bulk topological invariant, called spin-Chern number, which seemed to the contain 
more information about the edge structure. Later, however, it became clear that the 
spin-Chern number is unstable to deformations of the Hamiltonian system and that 
the only invariants for the spin-Hall effect are of Z2 type, rather than integer type 



Despite of the above concentrated efforts, that actually lead to a deep 
understanding of the topological insulators, there are still a few open problems. First, 
it is not completely clear how to describe the topological phases for non-crystaline 
systems. Secondly, it is well known [1] [H [31 [3] that the spin edge current is not 
quantized in the Quantum Spin-Hall systems and it remained an open problem to find 
an observable that has such a quantized edge current. It is worth mentioning that 
for the Integer Hall Effect, a complete, rigorous description of the edge states that 
goes beyond crystaline systems was achieved only in 2001 by Kellendonk, Richter and 
Schulz-Baldes [TJ] . The present paper was inspired by a later work of these authors [13] 
which deals with the quantization of edge currents for half-plane continuous magnetic 
operators in the presence of weak random potential. The formalism was put into an 
abstract setting in Ref. (TB], which actually provided the guiding lines for the present 
paper. This general formalism was applied in Ref. [IT to a simpler problem, namely 
the quantization of edge currents in Chern insulators with rough edges. The technical 
estimates derived in this paper are important for the present analysis. 

The present paper addresses both open issues mentioned above. Using the time- 
reversal invariance property of Spin-Hall systems, we define an observable and its 
corresponding current and we show that the expectation value (taken only over the 
spectrum in the insulating gap) of this current is quantized and that the quantization 
is described by the index of a Fredholm operator. This is our new topological invariant, 
which we call the edge index. For the model considered in Ref. [6] , we show that this 
invariant takes the same value as the Spin-Chern number defined in same reference. 
For the general case, we show that the edge index takes the generic values or 2, in 
line with the Z2 classification of the time-reversal invariant insulators. 

2. The model 

To be concrete, we consider non-interacting electrons on a honeycomb lattice (see Fig. 
1) described by the bulk Hamiltonian of Ref . [6 : 



|7l[9l[T0l[l2]. 
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Figure 1. The figure illustrates the honeycomb lattice, an example of paired 
sites with the index a attached to each site, and a random edge F. The contour T 
never crosses the bonds between the pairs and is contained between the vertical 
lines at — D and D. 
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Figure 2. The figure illustrates the energy spectrum of the Bloch Hamiltonians 
-fffe corresponding to the Bloch decomposition of the edge Hamiltonian relative to 
the translational symmetry along the homogenous edge. 

This particular model does not play any critical role in our analysis, except that it 
displays all the general features that we mention in the following. The Hamiltonian 
of Eq. [T| has time reversal symmetry and is a good model for electrons in graphene 
[2]. The first term is the usual nearest neighbor hopping term, the second term is an 
intrinsic SO coupling preserving the lattice symmetries and the third term stands for 
the Rashba SO coupling. For details about the notation please consult Ref. [F. In 
the following, we consider that we are in the Spin-Hall part of the phase diagram of 
the model [T]. 

The bulk model displays two top bands and two bottom bands separated by 
a gap. The two bottom bands have opposite Chern numbers c=±l, so their total 
Chern number is zero. When Vr=0, Sz commutes with the Hamiltonian and the 
model Eq. [l] reduces to a spin up and a spin down decoupled Haldane models [T5] . 
In contradistinction to the Chern numebr, the Spin- Chern number introduced in 
Ref. [6] is nontrivial: if Vr=0, it reduces to Cs = C| — (=±2 for the model Eq. fll 
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depending on the sign of Vso)- Cg can be generalized to the case when Sz is not 
conserved, like when the Rashba term is present. After extensive numerical analysis, 
Ref. [B] concluded that the Spin-Chern number remains quantized when Vr and a 
weak disorder are turned on. 

The special topological properties of the bulk energy bands have non-trivial 
consequences for the surface states spectrum when an edge is cut on a bulk sample. 
Let us briefly discuss the edge spectrum for a homogeneous edge. In this case we 
can use the Bloch decomposition with respect to the periodicity along the edge and 
write the edge Hamiltonian as a continuous direct sum of Bloch Hamiltonians Hk- As 
illustrated in Fig. 2, the spectrum of each Hk consists of upper and lower continuum 
parts plus two nondegenerate (excepting fc=0), discrete eigenvalues. These discrete 
eigenvalues for different fc's assemble themselves in two bands, shown in red color in 
Fig. 2. If the Rashba term is zero, one band corresponds to the spin up and the 
other band to the spin down. Thus, while the charge moves in opposite directions for 
these two bands (leading to zero charge current), the spins move in the same direction 
and consequently the edge carries a dissipationless spin current. The edge modes are 
protected by the time reversal symmetry, which means no gap can open in the edge 
spectrum, even when the Rashba term is turned on. While Sz is no longer conserved 
for this later case, the edge still carries a dissipationless spin current, thought no longer 
quantized. Because of the last fact, the theory of Quantum Spin-Hall is still missing a 
topological invariant that could tell how many edge bands one should expect in more 
complicated models. Finding such an invariant is the goal of the present paper. 

Our analysis will be done on an equivalent system, a triangular lattice with 4 
quantum states per site. This system is obtained by considering the honeycomb 
lattice as composed of pairs of sites sitting on a triangular lattice. For example, 
the 4 quantum states residing on the pair of sites circled in Fig. 1 can be thought as 
4 quantum states residing at a new lattice site positioned at the mid point between 
the pair. This way we obtain an equivalent triangular lattice model with 4 quantum 
states per site (see Fig. 2). The Hilbert space is now spanned by the states: 

|n, a),n = (^1,^2) G triangular lattice, a = (a, a), (2) 

where a=l,2 is the index introduced in Fig. 1 and a is the spin index. The triangular 
lattice sites are described by (ni, 712), where ni and 712 represent the coordinates along 
the two directions shown in Fig. 2. The bulk Hamiltonian becomes: 

^0 = EEFab'|n,a)(n',b| +fX'|n',b)(n,a|]. (3) 

n.n' a,b 

The coefficients r"jJJ can be computed from Eq. [l] but their explicit expression is not 
needed here. The flrst sum is over the nearest neighbors. 

We now consider the system with the edge. In the lab samples, the strongest 
irregularities are probably seen at the edges of the samples, so here we will concentrate 
at this type of disorder and we will neglect the bulk disorder (the mathematics 
still work for weak bulk disorder). On the honeycomb lattice we consider random 
contours F, like the one shown in Fig. 1, their main features being that they never 
cross the bond between the pairs and that they are confined within —D < ni < D, 
where D will be fixed from now on. On the triangular lattice, F can be described 
by a sequence {7n}n, where 7„ gives the deviation of F from the axis ni=0 at 
the row n2—n of the lattice, as illustrated in Fig. 2. We have 7„ G X, with 
J = {-D + 1/2, -D + 3/2, ...,D- 1/2}. Thus, F can be viewed as a point of 
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the set = T^°°: T—{. . . ,7_i,7o,7i, . . .}. On the set fi, we introduce the product 
probabihty measure, denoted by dT, which is the infinite product of the simplest 
probability measure on X: J f{n)du{n) — ^ SneJ /("■)' being any function 
defined on I. We remark that dV obtained in this way is ergodic relative to the discrete 
translations along the vertical direction of our lattice. We will use dV to average over 
all possible contours V. 

The system with the edge is defined on the Hilbert space Tir = 
span{\n, a), n to the right of F} and its Hamiltonian is given by Hp : Hr Ti-r, 

Hr^J2 EPaS' In, a) (n', b| + f |n', b) (n, a|], (4) 

n,n' a,b 

where the first sum is restricted to the sites located to the right of T. Hy remains 
time reversal invariant. 

3. The main result 

We define now the central observable. As it was pointed out for the case of 
translational, time reversal invariant, half-integer spin Hamiltonians [11 [H |S], the 
Hilbert space can be divided in two invariant subspaces. This remains true when 
the translational symmetry is broken. More precisely, the Hilbert spaces Tir can be 
decomposed as = TL-iT) ® 7Y+(r), where the orthogonal subspaces 'H±iJ') have 
the following special properties: 

6'H± (r) = (r) and Ht'H± {T)^n±{T), (5) 

where 9 denotes the time reversal operation, 9 — e^'^^y/''^K {K= complex conjugation). 
An important observation here is that the construction is not unique. Let us denote 
by Hj_ the orthogonal projectors onto 'H±{T) and define Sp = H!^ — HL, where we 
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reintroduced the index T to remind that the operator is defined on Tir- Our central 
observable is defined by the self-adjoint operator: 

Xr = i(y^S^ + S^y^), (6) 

where ?/r|n, a) = 7i2|n, a), defined on Tir, is the observable giving the vertical 
coordinate. The self-adjoint property of the central observable can be demonstrated 
by following a technique developed in Ref. [T^ . 

Our systems with edge and the observables Xr have very special properties under 
vertical translations of the lattice. Let 

Un\{ni,n2),a) = \{ni,n2 - n),a). (7) 

be the implementation of the lattice translations along the n2 direction. These 
translations can also be extended to a map t„ acting on the space of all possible 
contours T. The map i„ simply shifts a contour downwards by n sites. We now can 
list those special properties: 

(i) The family {7?r}rGn is covariant: UnHTU*^ — Ht„T- 

(ii) Based on 1, we can choose Sp such that ; = Ej^p. Moreover, [Sj,,i7r] = 0. 

(iii) The central observable obeys: 

u„Xr< = Xt^r + nS^r, [^r, ^r] = 0. (8) 

(iv) For any function /(e), commutators of the form [Xr,f{Hr)] form covariant 
families : 

u^[Xr, /(i?r)]< = [Xt,,r, /(fft„r)]. (9) 

We are now gearing towards the main result. We denote the spectral projector of 
Xr onto the spectrum inside interval [n— 1/2, n+1/2) by TTr{n). Note that, at least for 
a small Rashba term, the half-integer numbers are outside the eigenvalue spectrum 
of Xr- This can be shown via estimates on the resolvent of Xr using techniques 
developed in Ref. [19 . If troA = Tr{7rr(0)A7rr(0)}, we define the current of Xr as:|lB] 

Jr = tro {p(^r)^^|^| = itro {p{Hr)[Hr, Xr]} . (10) 

Here p(e) is the statistical distribution of the quantum states. Since we are interested 
in the contributions from the edge states, we assume that p(e) is a smooth function 
with support in the bulk insulating gap. 

Tight-binding Hamiltonians like were analyzed in Ref. [T7]- With the 
assumption that the amplitude of 7rr(0)|n,a) decays sufficiently fast for large \n2\, 
the technical estimates given in Ref. [17^ assure that, in the present article, all 
the operators appearing inside the traces are trace class (so the trace is finite and 
independent of the basis set used to compute it) and all the sums are absolutely 
convergent. 

Main Statement. Let F{e) = p(e). Note that F{e) is smooth and equal to 
1/0 below/above the bulk insulating gap; also F'{e)=—p{e). We define the following 
unitary operators: Ur = e~'^'^^^''^^\ If TTp is the projector onto the non-negative 
spectrum of Xy, then: 

dr Jr = ^ Ind {7r>;7r7r>} . (11) 
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This is our main statement. Let us comment on it first. The index is an integer 
number, defined on the class of Fredholm operators as: 

Ind^ = dim Ker[A] - dim Ker[A*]. (12) 

It has very special properties, the most important being the invariance to norm- 
continuous deformations of the operator that keep the operator inside the Fredholm 
class. In our case, it follows from the estimates of Ref. [T7] that, as long as the the gap 
remains opened and the support of p(e) remains inside the gap, we can deform p(e) or 
Hr without changing the index. Moreover, the index is independent of the contour F. 
To see this, we turn off the Rashba term (without changing the index) and reduce the 
system to two decoupled Chern insulators. But for Chern insulators, it was already 
shown in Ref. [17j that the index is independent of contour F. 

We now show that the index is equal to the Spin-Chern number introduced in 
Ref. [5]. We take F as a straight vertical line. Without changing the index, we can 
turn the Rashba term to zero. In this case the up and down spins decouple and we can 
take H± as the spin up and spin down invariant subspaces, respectively. Definitely 
Eq.[TT]applies equally well to the case when the set reduces to one point, the straight 
contour Fq (all we have to do is to take D—0). Then we have the following practical 
way of computing the index: 

Ind{^> C/ro7r> } = ztro{p(i7rJ[i/r„, XrJ} 

= ^Tr^{^T^MpiH^o)[H^„,y^oWM} (13) 
-zTrx {TTro (O)p(i?ro ) [Hvo , VToI^To (0)} 
Using the Bloch decomposition, this becomes 

TT 

E / iPi4.k)dkei,^p{ei,)dkeU]dk, (14) 

k——'n 

where e]j'^ are the edge energy bands. Since J p{£) — 1, each integral gives the 
difference between the number of forward and backward moving bands for the 
corresponding spin, known to equal the Chern number for the corresponding spin. 
Thus, the index is equal to the difference between the Chern numbers for spin up 
and spin down, i.e. it takes the same value as the Spin-Chern number introduced in 
Ref. IS]. 

Note that our main statement is about the average of the edge current and not the 
current itself. However, since the family {Hr}ren is covariant relative to translations, 
which act ergodically on f2, the spectrum of Hr is non-random. This implies that, if 
the edge spectrum becomes localized for a non-zero measure subset of fi, it will be 
locaUzed for all contours, except a possible zero measure subset of il. But this cannot 
happen, exactly because the average of the edge current is non-zero for Spin-Hall 
insulator. This allows us to conclude that the rough edge cannot destroy the edge 
conducting channels. 

4. Sketch of Proof 

With our assumption that the amplitude of 7rr(0)|n, a) decays sufficiently fast for large 
|n2|, it follows from the technical estimates of Ref. jlTj that ttj? J/rTTp is in the Fredholm 
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sign(x-3) - 
sign(x-2) — 
sign(x-l) 



sign(x) 

sign(x+l) 

— sign(x+2) 

- sign(x+3) 




Figure 4. A graphical representation of sign(x + n) = S(x). The top hnes 
represent the shifted sign functions sign(a; + n). The sum of the top lines results 
in the stair like function S{x) represented by the bottom line. 



class. Let TTp be the projector onto the negative spectrum of X-p and Sp = TTp — TTp . 
We compute the index using the formula: [121 HZ] 

Ind{7r>C/r7r>} 

= -5 E Tr{Mn)iU^ - /)Pr, C/r]^r(n)}, ^^^^ 

n 

where the sum is absolutely convergent. The projectors 7Tr{n) leave the subspaces 
'H±{T) invariant, so they decompose in a direct sum: 7rr(n) — TTp (n) ©TTp (n). Similarly 
for Ep: Sp = SjT © Sp . Due to property (3) listed above we have the following fact: 

u„7rp (m)M* = 7r^^p(rn =F "-)■ (16) 

We consider now the average over T. Since the index is independent of F, the operation 
can be omitted for the left hand side. On the right hand side, we use the fact that 
the trace of trace-class operators is invariant to unitary transformations and that the 
measure dT is invariant to the mappings t„, to write: 

-2Ind{7r>C/r7r>} 
= ^ / dr Tr{?/_„7rp (n)([/i^ - /)[Sp, Ur]7T^ {n)u*_J 

n 

+ ^/dr TrK7r+(n)(C/* - /)[Sp, [/r]^+(n)<} 

(17) 

(Tr{^-_p(0)((7;_P - /)[«_„Spul„, f/t_r]^r_„r(0)} 

+Tr{<p(0)(C/*^P - /)KS+<, ;7*„p]7r+p(0)}) 
= E /«fr(Tr{^p (0)({/i^ - /)[w_„E-^p7.1„, C/p]^p (0)} 

n 

+Tr{Ti+mU^ - /)KE+ ^^p<, ;7r]7r+(0)}). 
One important observation here is that: 

UinSf pM±„ = sign(X^ + n), (18) 
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(sign(a:)= the usual sign function) so we can draw the partial conclusion that: 

Ind{^>C/r^f?} = -i / dTti-oiiU^ - /)Esign(Xr + n), Ur]} (19) 

n 

As illustrated in Fig. 3, 

'sign{Xr + n) ^ S{Xr) (20) 



where S{x) is the staircase function shown in Fig. 3. But S{x) — 2x + s{x) where 
s{x) is a bounded periodic function s{x + n) = s(x). Based on this observation, we 
show that the contribution to the index from s(Xr) is zero. Indeed, we can follow 
Refs. [TBI [17] to show that, and under certain circumstances satisfied here, 

J dTtTo{ArBr}^ J dTtTo{BrAr}, (21) 

for any covariant operators Ar and Br leaving 7Y±(r) invariant. Since s{x) is bounded, 
we can open the commutator below, 

/drtro{(C/^-/)[s(^r),C/r]} 

= JdT troiiU^ - I)s{Xr)iUr - /)} (22) 

-JdT tro{{U^ - /)(C/r - /)s(Xr)} 



and s{Xr) is covariant since s{x) is periodic, so due to Eq. 21 the last two terms cancel 
each other identically. Thus, we arrived at the conclusion that: 

Ind{7r>C/r7r>} ^ ' j dV tro{(C/f^ - /)[Xr,C/r]}. (23) 

But this is exactly Eq. 42 of Ref. [17], with yr replaced by Xy- Thus we can repeat 
the steps of this work to complete our proof (note that property (4) is needed for this). 



5. Discussion 



Our construction is based on the splitting induced by the time-reversal operation 
9: 7ir = ^-(r) © '^+(r). This sphtting is in general not unique, but we introduced 
several constraints that limit the number of choices. These constraints are: the special 
requirement (ii) mentioned in Section 3 and the fact that the kernel of the operator 
Ep = n!|_— IIL needs to be rapidly decaying. This warrants that our assumption stated 
before the Main Statement holds true. Since the expectation value of the current is 
taken only over the states inside the bulk insulating gap, we have to consider only the 
splitting of the states inside this energy window. 

Let us first restrict our discussion to homogeneous edges, in which case the 
the system is mapped into itself by the discrete translations and the wavenumber 
k parallel to the edge is a good quantum number. Dropping the index F, which is 
no longer needed, the condition (ii) reads: WnSw* = S. Thus, the projections n!|_ 
are translational invariant, thus they must be given by sums over the k fibers. For 
example, the "+" projection must be of the form: 

nV = E (24) 

where 'ipn,k represent the Bloch functions. The summation goes only over a partial 
number of band indexes and the integral over k could in principle go only over parts 
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Figure 5. The diagram shows the bulk insulating gap and two surface bands 
crossing it. Both (a) and (b) situations show possible splittings of the states in 
this energy window into 'H± such that = Ti. — , but only situation (a) leads 

to projectors n!|_ = f dk |i/'±,fc} (V'±,fe I decaying exponentially with the separation 
l"2 - "al. 

of the Brillouin zone. Note that the phase of the Bloch functions are not relevant 
here, a good news because the phases are in general difficult to control. Now, since 
the kernel of this projector, n^(ni, 712; n[, n'2), must decay rapidly with the separation 
1^2 — n'2\, the integral over k must involve the whole interval [0, 27r] and the band "(An.fc 
must be analytic of k. For time reversal invariant spin 1/2 systems, the spectrum is 
at least doubly degenerate. For the case when the degeneracy is strictly two-fold, or 
quaternionically simple pOj, the bands in crystaline systems come in Kramers pairs 
and one can easily form the projectors P\. with the required properties. 

Let us exemplify. Consider first the model Eq. [T] As already discussed, there is 
a pair of Kramers bands crossing the bulk insulating gap. Fig. [5^ shows a properly 
chosen splitting, which gives projectors P^. that are exponentially decaying with the 
separation \n2 — n'2\- In contradistinction. Fig. |5|d shows a bad splitting, which leads 
to projectors decaying only as l/|n2 — njj. As already discussed, the unique choice 
shown in Fig. |5^ leads to an edge index equal to 2. 

We consider now a more complex situation in which we have more bands crossing 
the bulk insulating gap. Let us consider the situation of Fig. |6^. This is not 
quaternionically simple and we know that this case is unstable. The degeneracies 
at k—Q are protected by the time reversal symmetry, but the other two degeneracies 
will be split by small perturbations. The stable situation is shown in Fig. [6]3, which 
is quaternionically simple. In both cases there seems to be more than one possible 
valid splittings of the states. However, if we want to define the projectors PJ_ so that 
we go continuously (more precise analytically) from situation (a) to situation (b), the 
splitting can be done in only one way, by incorporating the bands that hybridize when 
the degeneracies are split into either or P^_. Thus, the only possibility of splitting 
the bands, for both (a) and (b) situations, is the one depicted in Fig. [6] Of course 
there is a freedom of choice in choosing the ± labels. With this unique choice, the 
edge index is zero. 

We can continue the argument for more complex situations, but we can already 
see the general conclusion: the edge index takes only the values or 2, in line with the 
Z2 topological classification of time-reversal invariant systems. Also, for homogeneous 
edges, the expectation value of the current of our observable is simply given by the 
charge current carried by the bands included in the -I- sector minus the charge current 
carried by the bands included in the - sector of the 6 splitting. This current is given by 
an index which is well defined if the 6 splitting leads to a kernel of that is rapidly 
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Figure 6. The diagram shows the bulk insulating gap and four surface bands 
crossing it. The case (a) is unstable and small perturbations lead to case (b), 
which is stable. The diagrams also show the unique splitting of the states in this 
energy window into 'H± such that 6^+ = Ti.-, splitting that gives projectors 
n!|_ = J dk |V'±,fe}(V'±,fcl decaying exponentially with the separation |n2 — rtjl 
and interpolating smoothly between situations (a) and (b). 

decaying with the separation \n2 — n'2\. So the message of our result is that, whenever 
such 9 sphtting exists and the index is non-trivial, there will be edge states that are 
topologically protected. 

The cases of a rough edge or when a weak bulk random potential is present 
are more complicated and, at this moment, we can only assume that the projectors 
n!|_ can be properly defined. There is already good progress in characterizing the 
edge states and computing the Z2 topological invariant for these cases [TU]. This work 
adopted an algorithm originally proposed by Fukui and Hatsugai [8] for computing the 
Z2 topological invariant for crystaline systems to the case of non-crystaline systems. 
We believe that we can adopt this new explicit computational algorithm to construct 
projectors n!|_ with the desired properties, for non-crystaline systems. 

6. Conclusions 

In conclusion, we found that the current of the observable X = ^[yS* -I- is 
quantized and that the quantization is given by the index of a Fredholm operator. For 
the model Eq. [l] this index was shown to take same value as the Spin-Chern number 
introduced in Ref. [6|. In general, the edge index takes the generic values and 2, in 
line with the Z2 topological classification of time-reversal invariant systems. 

Our result provides a non-trivial topological invariant that relates directly to 
the edge of the Quantum Spin-Hall system. The robustness of the edge modes 
to continuous, time reversal invariant deformations of the model can now be 
understood from the special properties of the index. We have made a fundamental 
assumption, namely that the amplitude of 7rr(0)|n, a) decays sufficiently fast for 
large |n2|. For homogeneous edges, we have shown explicitly how to construct S* 
with exponentially decaying kernels, in which case the fundamental assumption holds 
true. It seems reasonable to assume that one can complete a similar construction for 
non-homogeneous edges, in which case the analysis shows that the edge conducting 
channels are robust against random deformations of the edge. 
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